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Cosmic inflation driven by branes wrapping the extra dimensions involves Kaluza-Klein (KK) degrees 
of freedom in addition to the zero-mode position of the brane which plays the role of the inflaton. 
As the wrapped brane passes by localized sources or features along its inflationary trajectory in 
the extra dimensional space, the KK modes along the wrapped direction are excited and start to 
oscillate during inflation. We show that the oscillating KK modes induce parametric resonance for 
the curvature perturbations, generating sharp signals in the perturbation spectrum. The effective 
four dimensional picture is a theory where the inflaton couples to the heavy KK modes. The 
Nambu-Goto action of the brane sources couplings between the inflaton kinetic terms and the KK 
modes, which trigger significant resonant amplification of the curvature perturbations. We find that 
the strong resonant effects are localized to narrow wave number ranges, producing spikes in the 
perturbation spectrum. Investigation of such resonant signals opens up the possibility of probing 
the extra dimensional space through cosmological observations. 
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1 Introduction 



Cosmic inflation [1, 2, 3] not only sets the initial conditions for the Hot Big Bang cosmology, but 
also seeds structure formation in the universe by generating the primordial density perturbations. 
Conversely, the primordial perturbations are a powerful probe of inflationary cosmology. Upcoming 
precision cosmological measurements will further reveal details of the perturbation spectrum, which 
contains significant information about the underlying physics of the inflationary era. Since inflation is 
sensitive to Planck-scale physics, this provides a fascinating opportunity of experimentally studying 
physics at energy scales far beyond the reach of current terrestrial accelerators. 

In this paper, we explore the possibility that theories with extra dimensions can produce specific 
signals such as spikes and/or oscillations in the primordial perturbation spectrum. Infiation driven 
by branes wrapping the extra dimensions is investigated, and we will show that in the presence 
of the excited Kaluza-Klein (KK) modes (i.e. oscillation modes) of the brane along the wrapped 
directions, sharp signals in the perturbation spectrum are generated due to parametric resonance. 
The excitation of the KK modes can be triggered during infiation as the inflaton brane passes 
by localized sources or features along the wrapped inflationary trajectory of the extra dimensions. 
Since the mass of the KK modes are determined by the size of the wrapped cycles, for small extra 
dimensions the excited KK modes are heavy and thus oscillate, leaving imprints on the density 
perturbations. In this sense the wrapped inflaton brane scans the extra dimensions, whose output 
is, as we will show in this paper, sharp signals in the primordial perturbation spectrum. Moreover, 
the signals produced from the KK tower have a periodic nature in A;-space. Thus all these features of 
wrapped brane inflation offer us an opportunity of extracting information about the extra dimensions 
through cosmological observations. 

Our work is motivated by a class of models in string theory where inflation is driven by D- 
branes wrapped on cycles of the internal geometry, e.g., [4, 5, 6, 7]. We especially have in mind 
the construction of [6] where the tension of the wrapped D-brane sources the potential energy 
driving inflation, and the monodromy of the wrapped brane extends the fleld space to yield large- 
fleld inflation. However, let us stress that our work is not limited to such explicit construction, 
but discusses rather generic phenomenon that arise when inflation is driven by physically extended 
objects. 

From the effective four-dimensional point of view, the zero mode position of the brane plays the 
role of the inflaton, while the excited KK modes are heavy oscillating flelds coupled to the infla- 
ton. Perturbations with wave modes which resonate with the heavy fields' oscillation frequencies 
are enhanced/suppressed while inside the Hubble horizon, generating features in the resulting per- 
turbation spectrum. Effects of heavy fields during infiation have much in common with infiaton 
potentials with sharp and/or repeated structures, and have been the subject of extensive study, e.g., 
[8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26]. (See also [27, 28] for observational 
analyses of spiky modulations in the primordial spectrum.) Normally one finds that such effects are 
spread over a rather wide k range in the perturbation spectrum, since it takes a few e-foldings during 
inflation for the inflaton to come back to its attractor trajectory, or for the heavy fleld oscillations 
to damp away. However in this paper, we show that parametric resonance due to the oscillating KK 
modes can source sharp features that are localized to narrow wave number ranges of the density 
perturbations. This is due to the Nambu-Goto action of the brane (or the Dirac-Born-Infeld (DBI) 
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action for D-branes) yielding nontrivial couplings between the heavy KK modes and the inflaton, 
including those with the inflaton kinetic term. Such kinetic couplings will turn out to be extremely 
efficient in producing strongly resonant features. Furthermore, resonant amplification is highly sen- 
sitive to the oscillating amplitude of the KK modes. Since the excited KK mode oscillations are 
damped by the inflationary expansion, strong resonance happens only for a narrow wave number 
range, thus results in producing sharp spikes in the perturbation spectrum. In this sense, our work 
is not only a wrapped brane realization of the well-studied heavy field physics during inflation, but 
explores a new possibility for generating extremely sharp features in the primordial perturbation 
spectrum. 

This paper is organized as follows: First we explain the setup for the wrapped brane inflation 
we consider, and derive its effective four-dimensional action in Section 2. The dynamics of the 
brane, i.e. of the zero mode inflaton and KK modes, are discussed in Section 3. Then in Section 4, 
we derive the evolution equation of the inflaton field fluctuations, or equivalently, of the curvature 
perturbations. The resonant signals in the perturbation spectrum have different nature, depending 
on the amplitude of the excited KK modes. KK modes with small amplitudes give rise to oscillations 
with unique envelopes on the perturbation spectrum. This case, which we refer to as weak resonance, 
is analyzed mostly analytically in Section 5. For largely excited KK modes, the parametric resonance 
becomes extremely efficient and gives rise to spiky features in the perturbation spectrum. Such 
strong resonance is investigated in Section 6. Finally we present our conclusions in Section 7. In 
order to analyze effects from KK modes in the Nambu-Goto (or DBI) action for wrapped branes, in 
Appendix A, we derive the second order action for field fluctuations in general multi- field inflation 
with various forms of kinetic terms. Explicit forms of the functions that show up in the evolution 
equation for field fluctuations in wrapped brane inflation are laid out in Appendix B. In Appendix C, 
we give detailed analyses of the evolution equation for perturbations from weak parametric resonance. 

2 Effective Action of a Wrapped 4-Brane 

In this paper we would like to address effects of heavily oscillating KK modes on inflation driven 
by wrapped branes. For that purpose, let us consider the simple case where the inflaton brane is a 
4-brane that stretches along the external directions, while wrapped on an internal 1-cycle. Inflation 
is considered to happen while the 4-brane is moving along another direction of the internal manifold. 
We will suppose the six-dimensional part of the metric relevant to us to take the following diagonal 
form: 

ds^ = Qp^u dx'^dx'^ + g„ dr^ + gx\dX^. (2-1) 

Here = 0,1,2,3) are the external four-dimensional coordinates, r denotes the internal di- 

rection along which the inflaton brane moves, and A is another internal direction wrapped by the 
brane, which we consider to be compactified by A ~ A + 2tt. For an explicit realization of wrapped 
brane inflation, we have in mind the construction of [6] which considers wrapped D4-branes in ten- 
dimensional type IIA string theory. The detailed form of the internal metric will also be chosen 
based on [6] later on, however we should also stress that the phenomenon with KK modes that we 
will discuss are not restricted to the explicit string construction, but are rather generic features of 
wrapped brane inflation. (See also [29] which discusses wrapped brane inflation with light KK modes 
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along the wrapped directions.) 

We start by considering the Nambu-Goto action of the wrapped brane. Taking the brane coor- 
dinates {M = 0, • • • , 4) to coincide with x^^ and A, then the Nambu-Goto action of the wrapped 
4-brane with winding number p is^ 

//• /•27tp 
d^^^y-det{GMN) = -Ta J (fx J dA^- det(GMiv), (2.2) 

where T4 is the 4-brane tension and Gmn is the induced metric on the brane world- volume. Gmn 
is given by the metric and the brane position r(x^, A) as 

GuNdi^Ui^ = g^.^dxf'dx" + gxxdX^ + grr {d^^rd^r dx^'dx'' + 2d^rdxr dx^dX + {dxrfdX^] . (2.3) 

One can check that its determinant takes the form 

det(GAfAf) = g {g\\ + gxxgrrg^^d^rdur + grr{dxrf] , (2.4) 

where g = det(g(^i,). Thus far, we have been including derivatives of r up to all orders in the 
Nambu-Goto action. 

To make our discussion concrete, throughout this paper we suppose the Nambu-Goto action to 
be the main contribution to the inflationary Lagrangian, in other words, that inflation to be driven 
by the tension of the wrapped brane. We further assume the four-dimensional metric to depend only 
on the external coordinates, i.e. gfj,u = g^uix), and consider the internal metric to take the following 
form during inflation, 

grr = A^, gxx = 5^2, (2.5) 

where and are positive constants. Thus the basic picture is that the wrapped 4-brane starts 
off at some nonzero value of r, and drives cosmic inflation as it moves towards smaller \r\ due to the 
wrapped brane tension. (2.5) can be considered as an approximation of the internal metric during 
inflation, i.e., for the region of the internal space where the inflaton brane moves along. We should 
also note that this is equivalent to the case constructed in [6] where the extra six-dimensional space 
of type IIA string theory is compactified on nil manifolds [30]. (Though (2.5) is not exactly the 
geometry of the compactified nil manifold, monodromy of suitably wrapped D-branes extends the 
effective field range, i.e. allows large ranges for r, and realizes an effective metric of the form (2.5) 
in the large field limit.) The internal metric (2.5) allows large-field inflation for the canonically 
normalized zero mode position of the brane (p with the effective potential V oc (p'^^^, as we will soon 
see. 

In order to obtain the effective four-dimensional action, we now expand the 4-brane position in 
the internal space as 

00 

r{x'',X)= vnix") e'^'^/P (2.6) 

n=— 00 

where f„ = r_„,. In the following we use the abbreviations 

df-dg^ g^'^d^fd^g, {dff = g^^d^fdj, \df\^ = g^^d^fdj, (2.7) 

■^Though we also denote the brane coordinate along the wrapped direction by A as in (2.2), it should be noted that 
the brane coordinate A has periodicity A ~ A + 2np. 
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and we also assume 



r,ro>0, 1 + A'^idrof > 0. 



(2i 



Then, expanding the action (2.2) in terms of the (nonzero) KK modes, one obtains the effective 
four-dimensional Lagrangian {S = J d'^xC) as 



2'KpBTirQ — 

C = 

7 



1 + ^ 7'^'-(9ro • drn) + —l^\drn 



2 7' ^ n" ir 



,2 u |2 



2 S2p2 ^2 



7^ 4 

— {dro ■ drn){dro ■ dvn) \ + (cubic or higher in KK modes) 



, (2.9) 



where Yln^Q denotes the sum over all nonzero n, i.e. ^^"^=-00 n^O' ^'^'^ 7 defined as 



7 



(2.10) 



The wrapped brane inflation under consideration is basically a slow-roll one, and the inflaton brane 
does not go into the so-called DBI regime 7 3> 1 [31, 32], however it is important to keep the 7 
factor since it gives rise to kinetic couplings between the zero mode and KK modes, as can be seen 
in for e.g. the A^^'^\drn\'^ term. 

We have skipped terms including cubic or higher order in terms of the (nonzero) KK modes rn 
and/or their derivatives dvn- It can be checked that such higher order terms in the Lagrangian are 
smaller than the quadratic ones shown in the {} parentheses in (2.9) under the following conditions: 



\A\drof\, \A\dro-drn)\, \A\drnf 



A^ 

q2 p2 J.2 



«1, 



(2.11) 



where n in (2.11) represents any nonzero KK mode. Hereafter we restrict ourselves to this condition 
and consider interactions with the KK modes up to the quadratic order. Under (2.11), one may 
consider the term 7^A^(c?ro • drn){dro ■ dfn) in the second line of (2.9) to be negligible compared to 
n\ in the first line, however it should be noted that when focusing on their kinetic couplings to 
dro, they both are of the same strength. (Actually, we will see that these terms both give important 
contributions to the resulting KK-mode effects.)^ 
Let us now redefine the fields as follows, 



irpT^^A^B 



1/2 

2 D \ 3/2 



1/2 1/2 



Vq X 



Re(r„ 
Im(r„ 



(for n > 0) 
(for n < 0) 



(2.12) 



*To be precise, (2.11) is the sufficient condition for tlie quadratic terms of the most excited KK modes to be larger 
than any other cubic or higher KK terms. We also note that terms in the Lagrangian being small does not necessarily 
guarantee their effects on the equations of motion to be negligibly tiny. However, in this paper we simply drop the 
higher order KK terms based on the condition (2.11). 
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so that the action contains only real fields, and the zero mode (f) becomes canonical in the slow-roll 
limit in the absence of the KK modes. Then the Lagrangian is rewritten as 



1 V'n, 



2 ■ ^^i^<p■^^|J„ 
6 



+ ^ E I W • d^r.? + \§ - YJ^nd^ ■ di^n)] , (2.13) 



where again Yln^Q ~ l^^-oo n^o- have defined an effective potential 



V{(t>) = //i°/3</,2/3 with n 
and a dimensionless constant 



a. 



1 A'^ 



(2.14) 



(2.15) 



Note that = a?_„, and a„ denotes the positive root of (2.15). The 7 factor is now expressed as 



7=1 + 



V 



2\ -1/2 



(2.16) 



and the condition (2.11) for neglecting the cubic or higher order terms of the KK modes is trans- 
formed to 



{d<py 



V 



V 



V 



a. 



2< 
" ^2 



< 1. 



(2.17) 



In the absence of the KK modes, the Lagrangian (2.13) in the slow-roll limit reduces to the 
canonical form of Lj \J—g ~ —{d(f))'^/2 — V{(1)), which realizes large-field inflation. 



3 Dynamics of Wrapped Brane Inflation 

Now that we have the effective four-dimensional action (2.13), let us study the inflationary dynamics 
of the wrapped 4-brane. We start by discussing the homogeneous background. 

3.1 Homogeneous Background 

Upon discussing the homogeneous equations of motion, we fix the background metric to the fiat 
FRW: 

ds^ = -dt^ + a^{t)dx^. (3.1) 

Throughout this paper an overdot is used to denote derivatives in terms of the time t, and the 
Hubble parameter is defined hy H = a/a. Then the Einstein equation of the action (2.13) gives the 
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Priedmann equation 



37^(7^ -DE ^yS'^I- gl^n^n) , (3.2) 



as well as the evolution equation 



+ 7(37^ + 1)(7^ - 1) E (l^n + r8$'^n- I'^Mn] • (3.3) 



The full form of the equations of motion of the brane positions (/> and ipn are rather complicated, 
so let us partially write down the equations. The equation of motion of the zero mode (p is 

= j^ + 3H^+^-^^ + ..., (3.4) 

where • • • denotes terms containing the KK modes '0n and their time derivatives. This expression 
suffices for discussing the slow-roll dynamics of the inflaton cj). The interactions between (p and ipn 
will be analyzed in detail in the next section when we study perturbations. 

As for the KK mode ipn, the term —2Vja'^'ip'^/(j)'^ in the Lagrangian (2.13) sources the effective 
mass that is set by the length of the wrapped cycle, 

where we have ignored the 7 factor. Only writing down terms that are most relevant for us, the 
equation of motion of the KK mode ipn is 

= V^„ + 3//V'n + 4ya2|| + --- . (3.6) 
3.2 Slow- Roll and Heavy-Field Approximations 

For super-Planckian field values (j) > Mp, the action (2.13) can realize large field inflation with the 
potential V{(l)) (2.14) where (p 

is a nearly canonical inflaton fleld. If the KK modes ipn are excited 
during inflation, then given that their effective masses ~ Va'^/cp'^ are larger than H^, the KK modes 
would oscillate and leave resonant imprints on the curvature perturbation spectrum. 

We examine such case in this subsection and analyze the fields' dynamics. Specifically, we 
consider the case where the equations of motion (3.2), (3.4), and (3.6) are well approximated by, 
respectively, the slow-roll approximations 



3M^i?2 ~ V, (3.7) 
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and the heavy-field approximation 



3i?(/>~-^, (3.8) 



i^n^-^Val"^. (3.9) 



We further suppose that the order of magnitude of the ipn velocity is given by 

i^n-V^'^aJ^, (3.10) 

when averaged over the oscillation period. We should remark that throughout this paper, we use "~" 
to denote that the orders of magnitude of both sides of the equation are the same. For more precise 
approximations, we use "~" . Effects of KK modes on the curvature perturbations are discussed in 
Section 5 based on the above approximations, then in Section 6 we go beyond this case. 

Let us now show that under the approximations (3.7) - (3.10), the amplitudes of the following 
three parameters need to be sufficiently smaller than unity (We note that throughout the discussions 
we do not consider miraculous cancellations among terms in the equations.): 

Note that these are all positive parameters, since we are considering positive (p, cf. (2.8). 

Comparing with the full equations of motion, one can show that the sufficient condition for (3.7) 
and (3.9) to hold is 



< 1. (3.12) 



Here we remark that when we simply write, for e.g., r]n as in (3.12), then rjn represents the parameter 
for all nonzero n. Then, assuming e, Vn^n ^ 1, one can show from (3.3) that^ 

^ = 0{e) + 0{nnCn). (3.13) 

It should be noted that differentiating both sides of approximate relations does not necessarily 
provide similarly good approximations. Hence, in order to estimate the amplitude of (p, let us 
introduce a parameter 

3H<j) 



dV/d(t) 

From the equation of motion of (p, one can obtain 



1. (3.14) 



dV/d(i) 

Combining this with the time-derivative of (3.14), one arrives at 

k 



+ K = 0{e) + 0{nn). (3.15) 



^ -3K + 0(e) + 0(7?„), (3.16) 



'Expressions such as y — 0{x) are used to denote that \y\ < \x\ 
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which shows that |k| damps as the universe expands while \k\ ^ e, ry„. Hence one can conclude that 
the amplitude of k soon approaches 

At = 0(e) + 0(r/„), (3.17) 

and thus 

l^ = 0{e) + 0{vn). (3.18) 
Hep 

Therefore the term e^/^Cn4>Mp^^ fj,-^^/^ in (3.12) is estimated to be of size 0{e^Cn) + 0{er]nCn)- One 
can also show using (3.18) that the approximation (3.8) further requires r]n to be small. 
Moreover, since 

one can neglect the time variation of (/) in the approximation (3.9) and see that ipn at the leading 
order harmonically oscillates as cos{2V^^'^ant / 4>) ■ This validates the order of magnitude estimation 
of ipn in (3.10). (See also discussions around (3.26).) 

In summary, we have seen that when the approximations (3.7) - (3.10) hold, and given that there 
is no miraculous cancellation among the terms, then the following condition is satisfied, 

e, Vn, < 1. (3.20) 

This guarantees a„ ^ (l>/Mp 1. It should also be noted that the smallness of implies that the 
effective mass of the KK mode ipn is sufficiently larger than the Hubble parameter during inflation. 
For later convenience, here we lay out the order-of-magnitude estimations of various quantities: 

7-1, ^ = 0{e^) + Oier^n), = 0{e) + 0[r^nin), (3.21) 

^^^Oie)^Oi,.l (3.22) 

^ ~ ~ 1. (3.23) 

We also note that the condition (2.17) that was assumed in the previous section is rewritten in 
terms of the small parameters as 

(er/„e')'/', e, {eVninf' , r/nCn, r^nin <C 1, (3.24) 
where we have neglected the spatial derivatives. 

We end this section by writing down an approximate solution for ipn beyond the leading order 
equation (3.9), which will be used upon discussing curvature perturbations in the next section. One 
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can check that the ansatz of the form^ 



c \ 1 

^n = ^cosl2an^{t-U)+B\ (3.26) 

(here C, t^, and 6 are constants) satisfies the full equation of motion of -f/^n within errors of C(e, ^^), 
given that H{t — t-t) < 1, i.e., the time scale of interest is not much greater than the Hubble time. 
Here we note that upon estimating the size of the error, we have considered cosines and sines to be 
~ 1, and made use of (3.7), (3.21), and (3.22), but instead of (3.23) we have used the ansatz (3.26). 
In a similar fashion, one can also check that during H{t — t^) < 1, the ansatz (3.26) satisfies (3.10), 
and also ipn/Hipn ~ i^n/Hipn ~ l/'^n as in (3.23). Since the oscillations of the KK modes are quickly 
damped, the above approximate solution which is valid for a time range of ~ H^^ will be useful 
upon discussing effects from the KK modes. 



4 Effects of KK Modes on the Curvature Perturbations 

The KK modes of the wrapped brane can be excited during inflation when the brane passes by 
localized defects along the wrapped direction, such as other branes and "bumps" in the internal 
manifold. In this section we compute the curvature perturbations from wrapped brane inflation, 
considering the KK mode oscillations. We start by applying the general discussions in Appendix A 
to our action (2.13) and calculate the field fluctuations of the inflaton (p, which will be transformed 
into the curvature perturbations via the (5AA-formalism. 

4.1 Fluctuation Equation 

We are interested in the case where the nonzero KK modes ipn are heavy (i.e. m^j^ H^), hence 
we neglect their field fluctuations and only consider fluctuations of the zero mode inflaton, 

^ = 4>o + -, (4.1) 
a 

where (po is the homogeneous classical background. We work on flat spatial hypersurfaces, and 
hereafter we drop the subscript denoting the background. Then the second order action for the 

^This ansatz can easily be guessed from the next-to-leading order approximation of the equation of motion of ipn, 

V^n + 4Val ^ + 3//^„ - ^7/;„ - 8a^ ~ 0, (3.25) 

where the third term in the left hand side is suppressed by ~ ^„ compared to the leading (i.e. first and second) terms, 
while the forth and fifth terms are suppressed by ~ e. We note that (3.26) is a good ansatz independently of 5„ ^ e. 
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field fluctuations (A. 29) is (see Appendix A for the detailed derivation'') 

S2 = J drd^x !^-LAG^46'^diqdjq) + y (^M^^ - C^^ - 3^C^^^ + ^B^^ (^q' - | , 

(4.3) 

where r is the conformal time dt = adr, a prime denotes a r-derivative, and the quantities except 
for q are those of the homogeneous background. The definition of the functions -B^^, C,^,/,, etc. are 
laid out in Appendix B. The second order action gives a linear equation of motion for q, which, after 
Fourier expansion 

q{r,x) = J^I d^ke'^^-^quir), (4.4) 

takes the form (cf. (A. 35)) 

= + l^agL + l^^^fc^ - - + zM^I±^a' + Z^^:!L±^a\ q,. (4.5) 
B(f,^ I 1500 a i500 ii00 J 

Here, k = \k\. We remark that is positive throughout the cases discussed in this paper. 
4.2 Analytic Study of Field Fluctuations 

In the rest of this section we carry out analytic calculations of the field fluctuations by restricting 
ourselves to the case studied in Section 3.2, where the excited KK mode amplitudes are small (cf. 
(3.20)) and the background field dynamics are approximated by (3.7) - (3.10). We will compute the 
leading order effects from the KK modes on the curvature perturbation spectrum. 

4.2.1 Evolution Equation for Small KK Excitations 

Under the condition (3.20), the homogeneous functions in (4.5) (cf. (B.4) - (B.7)) can be evaluated 



as 



B,, = l-2^al§ + l^^ + Oie), (4.6) 
^ = 0{e') + 0{er,^^r.), ^ = 0{e') + 0{er^r.), ^ = 0{e), (4.9) 



^Discussions in Appendix A are applied to our case by substituting 
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where we have used (3.21), (3.22), and (3.23) upon estimating the amplitude of the dropped terms. 
For the term Bff,^/H (4.7), the first two terms in the right hand side each has amphtude ~ 7]n. 
We also note that for i?^^ and L^G^^, the second and third terms on the right hand sides are of 
order r/„^„, and the other dropped KK terms (terms depending on tpn and "0™) in (4.6) and (4.8) are 

The conformal time can be computed by integrating 



d 



1 



1 + 



H 
IP 



dr. 



Hence one finds 



1 



(l + 0(e) + 0(77„e„)). 



where we have took r to (at the leading order) approach from the negative side as a 
one can further show 



a 
a 



1 



(1 + 0(e) + 0(7?nCn)) 



a 
a 



(l + 0(e) + 0(r/„en)). 



(4.10) 



(4.11) 



oo.^ Then 



(4.12) 



Thus by combining the above estimations, one can rewrite (4.5), keeping only the leading order 
contributions from the KK modes as 



= 4 + 



VH 



Ik 
T 



+ 



1 



2 + -Y.'^l 



VH 



^. (4.13) 



Let us repeat that the ipn terms explicitly written inside the parentheses on the first and third lines 
have amplitude of ~ r/„, while that in the second line is ~ ?7n^n-^ One may except the term ipn/^ in 
the second line to have much smaller contribution than the other KK terms, however this is not the 
case for the following reason: The field fluctuation experiences parametric resonance with the KK 
mode oscillations when their frequencies synchronize. This happens when k/a ^ Jn-KK ~ OinV^^'^ 
i.e., when k V2 ~ £-2. 

(We will see this explicitly in Section 5.) Thus for that experience 
parametric resonance, the term k'^r'^tp'^/V is of order rjn/£,n until the wave mode passes the resonance 
band. Thus in the end, this term gives contributions comparable to those from the KK terms of 
order r/„ in the first line of (4.13).^*^ We also note that the KK terms in the third line turn out to 
be irrelevant inside the horizon, as we will soon see. 



Strictly speaking, the 0{e) and 0(77,1^,1) in (4.11) are integrated values of the parameters over a finite period of r, 
however for simplicity we treat them as 0[e, rinS,n)- 

^We also note that the term i}^/V in the second line of (4.13) is the leading KK mode contribution to the coefficient 
of fc^T^, since the next-to-leading KK terms are of ©(er;™^™, r;ij^jj). 

^''One may then naively think that the term k^r'^ip^/V should have much larger effects than the other KK terms, 
however their contributions turn out to be comparable when taking into account that it is actually the absolute value 
of Qk that matters for physical observables. 



12 



We stress that the leading KK terms shown expHcitly in (4.13) (or in (4.6), (4.7), and (4.8)) 
arise from the foUowing three terms in the Lagrangian (2.13) through their kinetic couplings with 
the inflaton, 

-27^E«'S' "iE^^v^")'' ^Ew-^^«)'- (4.14) 

n^O n^O n=/=0 

These kinetic coupling terms play important roles even when we go beyond the small KK excitations 
and study strong resonance in Section 6. 

We expand the field fluctuation in terms of rjn as 

Qk = q'j^^+q^^ + --- , (4.15) 

where = 0{r]'!^), and solve (4.13) at each order. (We remark that here we are implicitly assuming 
a hierarchy between e and rjn (i.e., either e <^ ijn or e ^ r/„), otherwise we may have to analyze the 
0{e) correction at each order in r]n-) Then the r/^ order (4.13) is, at zeroth order of e, 

C+(*=-^)«r-0. (4.16) 

The r/^ order (4.13) is, at zeroth order of e and 

where in the right hand side we have used (3.9). 
4.2.2 Zeroth Order in rjn 

Let us first focus on the rj^-order equation (4.16). The general solution to this equation is given by 
a linear combination of the Hankel function 

^<(-M = /S(^-l)^-'' ("8) 

and its complex conjugate. The explicit form of q^^ is determined through quantizing the field 
fluctuations. We stress here that the KK modes are assumed to be excited during inflation, hence 
in the early stage of inflation their effects are absent, i.e. qk = q^/^^ ■ 

From the action S2 = J dTd?xC2 (4.3), one obtains the conjugate momentum of q, 



We expand q in terms of the mode functions q^ and further assign annihilation and creation operators 

4 



(ofc and at respectively). 



Qke ^'"■^fflfc + (qke 



13 



and impose the following commutation relations 

[afc, at] = (2^)35(3) 



[fflfc, Oq] = K, ay = 0. 

as well as 

[g(r, x), q{T, y)] = [n(r, x), n(r, 2/)] = 0. 

Requiring to realize the positive frequency mode of the plane wave solution oc e^**^"^ in the 
early times when the mode is well inside the horizon, and taking into account that q^ is initially 
equivalent to q^^\ then we can choose the solution (4.18) as the mode function, i.e. qk = q^^^ = 
CyJ —TH^)j—kT), at times before the KK modes are excited. We choose the normalization C (which 



'3/2^ 

is a constant at zeroth order in e), such that the commutation relations (4.21) and (4.22) are satisfied. 
Hence C can be taken as 

where the estimated error 0(e) arise also from B^^ (4.6) and a' /a (4.12). (The error 0{rin£,n) is 
absent here since we are discussing times before the KK excitations.) 

After the KK modes are excited, the solution (4.18) of qu is succeeded to q^^\ and KK mode 
contributions q^j^\ q^^\ ■ ■ ■ arise. Here, let us note that when starting with an initial condition that 
is independent of the direction of k (such as (4.18)) and satisfies the commutation relations (4.22), 
then one can show that (4.22) is satisfied as long as the mode function qj. follows the equation of 
motion (4.5). 

Therefore, within errors of 0{e), we have obtained q^^ as 

M 1 1 f ^ ,\ ^-ikr 



4.2.3 First Order in ri„ 



Let us express the source terms in the right hand side of (4.17) as functions of the conformal time. 
In order to analyze the leading effects from the KK modes, it basically suffices to obtain the leading 
order expressions. By solving (3.7) and (3.8) using r ~ —1/aH, one finds 

,/)2:^02 + ^M2ln-, (4.25) 

where the subscript * denotes values at some fixed time r*. For the KK mode we use the expression 
(3.26), giving 

~ /a.\3/2 f yi/2 1 
V^^-V'n*^— j cosl2an — {t-t,)+en>, (4.26) 



^^Strictly speaking, the error 0(e) in (4.23) is actually e further multiplied by factors such as k and r. However, 
here for simplicity we ignore such factors. 
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where ipn*, a*, t*, and On are constants. We obtain the expression for ipn by differentiating the 
ansatz (4.26), which is at the leading order (for H{t — t*) < 1), 



i^n ~ -'^an——'4^n* y — j siu i2an——{t - + On 



(4.27) 



Moreover, given that |-^ln(^)| <C 1, one can check that 



dt d <^ —— ln(— a*r//) 



Integrating (4.28) and further using (3.7), we obtain 



(4.28) 



/i5/3^1/3 r. 



1/3' 



(4.29) 



Hence by combining the results, one can write the source terms of (4.17) as functions of the 
conformal time, 

16 2 VWVi 



~ -16 • 3 



1/2 ' T 



E 



M„ . / T 



1/3N 



E 



l_cos -4.3"V,^l.i(- 



1/3N 



+ 20„ 



+ 20„J>, (4.30) 
(4.31) 



Here, in order to avoid clutter we have left 0, which is given in terms of r in (4.25). 



Introducing the ratio 



~ (0)' 

9fe 



and substituting (4.24) for the equation (4.17) can be recast in the form 

"(1) , nf ^ , 1 , 7(1)-, ik 16 2 '^ntp: 



(4.32) 



(4.33) 



.(1) 



Supposing, for instance, that the KK modes are suddenly excited at time Texc, then one can obtain s), 
by solving (4.33) (with (4.25), (4.30), and (4.31)), with the initial conditions s^^''(Texc) = s^^\tcxc) = 0. 
We choose the vacuum as 

afc|0)=0, ^fc, (4.34) 



^^Here we use the approximate solution (4.26) beyond the leading order equation (3.9), since the damping of the 
oscillations tpn oc a~^^^ need to be taken into account for understanding the behavior of the KK effects. 
^■^Here we suppose that the integrated errors stay small, as in Footnote 8. 
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hence the two point function of the field fluctuations is 

(2vr) 



{0\Scl,kS^^,\0) = + k') Iq^f . (4.35) 



a 



Considering the KK mode effects to first order in r/„, then 

P {l + 2Re(^))} - (^^^ + l) {l + 2Re(4^))} . (4.36) 

4.3 Curvature Perturbations 

Finally we use the JAA-formalism [33, 34, 35, 36] to compute the curvature perturbation which is 
expressed as 

dJ\f 

Ck ^ -Q^Hk (4.37) 

at the leading order of the inflaton field fluctuation. Here, N is the e-folding number between an 
initial flat hypersurface and a final uniform density hypersurface. Let us repeat that, even though 
there are multiple degrees of freedom involved in wrapped brane inflation, the KK modes are massive 
hence the curvature perturbations are sourced by the inflaton field fluctuations. We take the final 
uniform density surface at some later time when the wave mode is well outside the horizon so the 
separate universe assumption is a good approximation, and also when the KK modes have sufficiently 
damped away such that the infiationary universe can be considered as nearly single-component. 
Computing (4.37) at such time, one obtains 

— ~ -^5</.fc, (4.38) 

0(p (f) 

up to corrections from the damped KK modes. Here H/cp can be considered as a function of cj) since 
now infiation is nearly single-field. 

The (damped) KK corrections, i.e. corrections to the single-component treatment, can be esti- 
mated from the Friedmann equation (3.2) which can be written as 

3M2f2 = V{1 + 0{e) + 0(r?„en)} , (4.39) 

and the equation of motion of (/>, 

3F(/) = -y'{l + 0(e) + 0(??„)}. (4.40) 

From these equations, the KK corrections to the expression (4.38) can be estimated to be of size ~ rj^ 
at the time when evaluating its right hand side. (Here, note that r/n decreases in time as the KK 
oscillations are damped.) On the other hand, the KK mode effects are also imprinted on the inflaton 
fleld fluctuations 6(j)k- We will see in the following sections that these effects on 5(j)k are at least of 
order ~ (??nCn) |exc which are values at the KK excitation. Therefore, by evaluating (4.38) at times 
when the wave mode of interest has exited the horizon, i.e. 

- A:t < 1, (4.41) 
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and also when the KK modes are damped such that the parameter rjn is smaU enough to satisfy 

Vn < {Vn^n)Uc, (4.42) 

then one can safely ignore the KK modes at around the final uniform density hypersurface but still 
capture the main KK effects through the inflaton field fluctuations. 

Combining (4.35) and (4.37) - (4.40), then at zeroth order of e and r/„ (we stress that this r]n 
is the value after the KK modes have damped away), one obtains the two point function of the 
curvature perturbations up to linear order in the inflaton field fluctuation as 

iCtCk') ^ {27rf + (^J^y (4.43) 

where the right hand side should be evaluated at times when (4.41) and (4.42) are satisfied. Further 
using (4.11)^^ and the order rjnexc expression (4.36), the two point function can also be expressed as 

iaC,') - {2nfs('Hk + fc')^ X ^^J^,y,, (1 + k'r') {l + 2Re(.«)} . (4.44) 

Here we should remark that the expression (4.44) does not capture the large (super-Planckian) 
variation of during inflation, since in the previous subsection we have focused on the KK mode 
effects and dropped 0{e) corrections upon calculating qk- 
Defining the power spectrum V(^{k) as 

r2 
fc3 



iCkCk') = i27Tfd^'\k + k')'^V^{k), (4.45) 



one sees that effects from the excited KK modes are represented by Re(s^^^) evaluated at a fixed 
time r when the wave number range k ^ k + Ak of interest satisfies (4.41), i.e., 

^^l + 2Re(4^)). (4.46) 

Here, V^oik) denotes the curvature perturbation power spectrum in the absence of the KK modes. 



5 Weak Resonance from Small KK Excitations 

We now analyze KK mode effects on the curvature perturbation spectrum by solving the evolution 
equation derived above. In this section we consider cases where the condition (3.20) holds and the 
inflationary dynamics are well approximated by (3.7) - (3.10). Such small KK excitations give rise 
to weak parametric resonance, sourcing oscillations on the perturbation spectrum. In this small KK 
regime we can carry out analyses analytically, enabling us to understand the particular oscillatory 
forms the KK modes source on the power spectrum, cf. Figure 2. We will go beyond the small KK 
approximations in Section 6 where we find strong resonance sourcing sharp spikes, but the basic 
properties of the KK signals will still be described by the analyses given in this section. 

^''As we have stated in Footnote 8, the estimated errors of r in (4.11) are actually those parameters integrated over 
a finite time period. However, since we fix r — >■ 0~ as a — >■ cxa, the integration is from r to the future and thus errors 
in r at the final hypersurface are not sourced by KK modes at times when they were excited. 
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5.1 Solving the Evolution Equation 

Recall from the discussions in Section 4.3 that for small KK excitations (3.20), their leading order 
effects on the curvature perturbation spectrum is obtained by computing the asymptotic value of 
2Ke{s^i^^). Focusing on wave modes k that are well inside the horizon when the KK modes are 
excited, then the evolution equation for s^^^ (4.33) prior to horizon exit (i.e. during —kr S> 1) is 
approximated by 




where (j) is given in (4.25). (Here one sees that the KK terms in the third line of (4.13) is irrelevant 
when inside the horizon.) The homogeneous solution of this equation is 

s^lir) = Co + C^e'^'\ (5.2) 

denoting oscillations with frequency k/ir and an offset. Since (5.1) is a linear equation, let us look 
into s^^^ induced by each of the source terms in its right hand side, and then add up the results at 
the end. We suppose that the KK modes are excited at time Texc, and set the initial conditions as 

4'^(w) = 4'^(rexc) = 0. (5.3) 

Hereafter we use the subscript "exc" to denote values at the KK excitation. 



5.1.1 Non-Oscillatory Source 



Upon dealing with the non-oscillatory source term (i.e. the one without the sine in the right hand 
side of (5.1)), we ignore the logarithmic r-dependence of (p (cf. (4.25)) and solve the equation 



-2r 



T 



/)2 
exc 



(5.4) 



where we have taken the arbitrary time r^, in (5.1) as at the KK excitation. Solving (5.4) with the 

initial conditions (5.3), then the real part of s^^^ which is relevant for the curvature perturbations is 

obtained as _ 

,9 r / 

r 



^exc 



COs{2A;(t - Texc)} 



(5.5) 



This solution induced by the non-oscillatory source consists of an oscillatory part with frequency fc/vr, 
and an decaying offset. Since we are now focusing on wave modes inside the horizon, in the above 
solution we have dropped oscillatory terms and (decaying) offset that are suppressed by Od/crl^^) 
compared to the shown terms. 
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5.1.2 Oscillatory Source 

In order to analyze s^^^ induced by the oscillatory sources, we focus on time intervals of order the 
KK-mode oscillations (which is much shorter than the Hubble time due to ^ !)• Then we can 
ignore the time dependence of the source terms except for that showing up explicitly in the sines, 
thus obtain 



-4«(r) + m4\r) ^ 16 . sin 



-4.3VV^ln(-l+20, 



(5.6) 



16 • 3^/^— ^ sin ((i^„^,r + @n*) 



^ sin (uJn*T + Qn* + ^ 



Here, the parameters with the subscript * represent their typical values during the time scale of 
interest at around some arbitrary time t = t^, and especially ipn* denotes the oscillation amplitude 
of the KK-mode ipn- In the second line we have expanded the log term and introduced the parameters 



(5.7) 



These parameters can be considered as constants during the time intervals of the KK-mode oscilla- 
tions, however at larger time scales (i.e. Hubble time or more) they vary approximately as 



rjn oc~ T 



(5.^ 



From the homogeneous solution (5.2), one can expect parametric resonance to happen when the 
oscillation frequency Un/^ii becomes equal to the characteristic frequency of the mode /c/vr, i.e. 
2k ~ a;„. This is at a time 

- Tres{k) = —r- > 7, (5.9) 



k^n k 



when the mode is well inside the Hubble horizon. 



We collectively study the oscillatory terms in (5.6) by analyzing the following equation 



4^\t) + 2iksl''{T) = iAsm{u:T + 9), 



'(1) 



(5.10) 



where ^ is a complex constant, while uj and are real constants. The particular solution of (5.10) 
is 

A ( Oh ^ 

(5.11) 



UJ 



A (2k 

< — cos(a;r + ©) + z sin(a;r + 0) 



2 - 4A:2 [ UJ 



for / Ak'^, and 



^i(-) = ^ {2.fere^(^^^+«) - cos(2fcr + 9)} 



(5.12) 
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for Lo = 2k. The particular solutions (5.11) and (5.12) connect in the a; — )• 2A; limit up to the 
homogeneous solution (5.2): 

While the homogeneous solution (5.2) sets the offset as well as oscillations with frequency k/ir, 
the particular solution (5.11) gives oscillations with frequency u/ln. Here, the parameters such as 
A and actually vary in cosmological time scales in the original equation for s^^\ However, since 
their varying time scales are much longer than the oscillation period ~ uj~^ (and also than k~^), 
one can consider (5.11) to be the solution at each cosmological time scale. In other words, the 
particular solution denotes oscillations with time dependent amplitude and frequency, in contrast to 
the homogeneous solution giving oscillations with constant amplitude and frequency. 

As for the resonant solution (5.12), the first term in the right hand side represents parametric 
resonance. The degree of parametric amplification is determined by how long the wave mode stays 
in the resonance band oj ~ 2k, see discussions in Appendix C.2. 

Detailed solutions for the equation (5.10) are given in Appendix C (e.g. (C.IO), (C.ll), and 
(C.12) for time dependent parameters), but here let us give a rough description on how the real part 
of s^^'' reacts to the oscillating sources of (5.6). 

Upon the KK excitation, if the KK-modes' oscillation frequencies are larger than the character- 
istic frequency of the wave mode, i.e. 2k <C Wnexc, then Re(s^^'') is insensitive to the KK oscillations 
and is well described by the homogeneous solution (5.2). The oscillation amplitude and offset are 
set by the initial "kick" at the KK excitation, and are of order ?/nexc?nexc- Such wave modes do not 
experience parametric resonance and eventually exit the horizon. 

On the other hand, wave modes satisfying 2k ^ u)ncxc go through the resonance band before 
exiting the horizon. When the KK modes are excited, oscillations of Re(s^^^) with frequencies k/ir 
and a;„/27r both can be generated with amplitude of order rjnexc^nexc- However the former keeps a 
constant amplitude (corresponding to the homogeneous solution (5.2)), while the latter decays as 

as the KK mode oscillations are damped (corresponding to the particular solution (5.11)). The 
offset is also of order r/„excCnexc- 

The wave mode undergoes parametric resonance with the KK oscillations while 2/c ~ a;„, during 
which the oscillation amplitude of Re(s^^^) is amplified. At this time Re(s^^'') oscillates only with 
frequency k/ir. The wave mode stays in the resonance band for Ar ~ — Cnexc7"res (cf- (5-16)), and by 
the time it leaves the resonance band, the oscillation amplitude of Re(s^^^) is parametrically amplified 
to 

~ fJnexcCnexcikn/ k)'^ . The amplification is peaked at a wave number ~ kn experiencing resonance 
right when the KK modes are excited. For larger k, the amplification is less significant since such 
wave modes enter the resonance band at later times when the KK oscillations have damped away. 
We also note that the fc/vr oscillation of Re(s[,^'*) with constant amplitude rjnexcCnexc (denoted by the 
homogeneous solution) survives through the resonance band. This can be the dominant oscillations 
for large k modes that are not so amplified at the resonance band. 

After leaving the resonance band (i.e. when 2k <C w„), then Re(s^^^) becomes insensitive to the 
source terms and maintains its oscillation with frequency k/n and amplitude at the resonance band, 
until the wave mode exits the horizon. 
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5.1.3 After Horizon Exit 

The evolution of s^.^'' after horizon exit can be analyzed in a similar fashion, by studying the equa- 
tion (4.33) for —kr <^ 1. (In such case, the second source term ipn/V in (4.33) is much smaller than 
the first term and thus can be ignored.) Analyzing the super-horizon version of (5.6) by neglecting 
the time-dependence of the parameters such as a;„, then one can easily check that the solution out- 
side the horizon consists of a constant term and decaying terms. Hence the oscillation of Re(s^^^) 
becomes frozen as the wave mode crosses the horizon, i.e. at —kr ~ 1. One can also check that the 
induced constant term for wave modes that are super-horizon at the KK excitation is tiny compared 
to the typical amplitude of Re(s^^^) for wave modes inside the horizon at Texc- 

5.2 Approximate Solutions 

Now let us present approximate solutions to the evolution equation (4.33) with initial conditions 
(5.3), based on the discussions above and in Appendix C. The results will be written in terms of the 
following parameters: 



'/nexc — /•^ 



2 • 3^/2 

kn^ , (5.14) 

Texcsn exc 



where V'n is the oscillation amplitude of the KK mode ipn, and kn denotes the wave number in the 
resonance band when the KK modes are excited, i.e. 2kn = ujnexc- Furthermore, since a;„ cx r"^, 
the time when the wave number k enters the resonance band can be estimated as 

W.)^<^^-|^. (5.15) 

"'Snexc 

The time range a certain wave mode stays in the resonance band is derived in (C.8), which can now 
be rewritten as 

Tf{k) - rres(fe) -eJLr,,,{k). (5.16) 



The approximate solutions for Re(s^ ) are obtained by adding up the contributions from the 
oscillatory sources (i.e. (C.IO), (C.ll), and (C.12) for (5.6)) and non-oscillatory sources (i.e. (5.5)). 
Here we show fluctuations induced by a single KK mode n, but the full expression is simply a sum 
of the following expressions over all excited KK modes. 

For wave modes that are sub-horizon at KK excitation but do not cross the resonance band, i.e. 
— 1/rexc < k < kn, one finds 

"^^(■^fc (^)) ~ ^nexc^nexc 

for Texc <T < - J. (5.17) 



2cos26'„ - M cos{2A;(r - Texc)} - 2cos26'„ + ^ ( — — 

/ \ ^exc 
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For modes that do cross the resonance band, i.e. kn < k, the solution before crossing the 
resonance band is: 



3 r 4.31/2^ 

cos < 1 m 



n cxc 



+ 29^ 



~\~ fjn exc 



1 1 / T 

- (cos 26n - 1) COS {2k{T - Texc)} - 2 COS 26*^ + - I 



for Texc < T < Tresik), (5.18) 



and after crossing the resonance band: 



Re(4'^(r)) ^ 33/2^,exc e^exc ( y 1 COS <! 2A;r + 



kr. 



4 • 

■^n exc 



1 + ln 



/n exc exc 



- (cos 29n - 1) COS {2k{T - Texc)} " 2 COS 26'^ + 
Z Z \ Tp 




1 



for rres(/c) <r<-j, (5.19) 

where for simphcity we have supposed that the parametric amphfication happens suddenly at r = 
Tj-es{k). Upon obtaining the first line of (5.19), we have used T^res — ??exc(^n/fc)^) which follows from 
fjn oc~ r^, (5.14), and (5.15). 

In the above results we have dropped sub-leading contributions to oscillations with frequency k/ir 
that would follow from simply adding the results (C.IO), (C.ll), and (C.12), induced by the two 
oscillatory source terms in (5.6).^^ We should also remark that since we have collected leading 
contributions from each source term, the above results can become inaccurate when the main terms 
in the expressions exactly cancel each other (for e.g., when 2cos20„ = 1/2 in (5.17)). 

When the wave mode exits the horizon at around r = — l/Zc, then Re(s^^'*) freezes out for both 
(5.17) and (5.19). 



In Figure 1 we plot the time evolution of Re(s[,^^), comparing the approximate solutions with 
numerically computed results. Here the inflationary parameters are chosen as i?i>exc ~ 8.2Mp and 
fj, ~ O.OOlGMp, such that the KK modes are excited at about 50 e-foldings before the end of inflation, 
and the COBE normalization value ~ 2.4 x 10^^ [37] is realized at the wave number kp which 
exits the horizon at the KK excitation. We have assumed the KK mode n = 1 with ai = 100 to be 
excited as ipnexc = O.OOlMp and 6i = (i.e. ipnexc ~ 0). This set of parameters is translated into 
the small parameters ^lexc ~ 0.082 and f/iexc ~ 0.0018, and realizes the resonance peak at around 
ki ~ 42/cp. We show two plots for k ~ 14/cp, 140A;p, each representing wave modes that do not/do 
cross the resonance band, respectively. The blue solid lines are obtained by numerically solving 
(4.33) (with (4.25), (4.30), and (4.31)), while the red dashed lines show the analytic estimations 
(5.17) - (5.19). 



'oxc?7ncxc^LxcSin26l„sin{2A:(r - Toxc)} in (5.17), and -4 



^^In particular, the dropped terms are —4:^^3~^^^kTc: 

3^^^k~^T^lfinaiccSm29„sm{2k{T — Texc)} ln (5.18) and (5.19). They can become comparable to the leading kf-K 

oscillatory terms for k kn, however, for such wave modes the parametrically amplified term (i.e. the first term in 
the right hand side of (5.19)) is basically dominant anyway. 
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Re(4' 




Figure 1: Time evolution of Re(s^^'') in terms of normalized time kpT. The KK mode is excited at 
kpT ~ — 1. Blue solid lines show numerically computed results, while the red dashed lines are the 
analytic estimations (5.17) - (5.19). Right (Left): Case for a wave mode that does (not) cross the 
resonance band. 



One sees that the analytic estimations well describe the numerical results, except for at around the 
time Tj-es in the right figure, where we have made the simplifying assumption that the amplification 
happens instantaneously. It should be noted that parametric amplification actually happens roughly 
as OC (r - Tres), cf. (C.6). 

We also note that in the analytic estimations, we have simply frozen the fluctuations at r = —1/k, 
which also give rise to errors for estimating the asymptotic values Re(s[.^^)|^_j.o- • However, such 
errors only lead to shifts in the phase of the oscillations in the Re(,s^^^) spectrum in /c-space, as we 
will see in the next subsection. 



5.3 Curvature Perturbation Spectrum 



The effects from the KK-mode oscillations on the curvature perturbation spectrum are captured 
by the asymptotic value 2Re(s^^'')|T-_!.o- , cf. (4.46). We obtain this analytically by freezing the 
expressions (5.17) and (5.19) at horizon exit r = —1/k. 
For k < -1/Texc: 

.0--0. 



Re(4^)) 



For -1/Texc <k<kn: 



For kn < k: 



2 COS 2^„ - M COS {2kTexc + 2) - 2 cos 29n - W7T— 



/2 

cxc , ^ 



krj 



COS 



4 • 3^/^ 

cxc 



1 + ln 



kn 



29n 



~l~ fjn exc exc 



- (cos 29n — 1) COS {2kTcxc + 2) — 2 COS 26„ 



(5.20) 
(5.21) 



(5.22) 



Here in (5.22) we have dropped a contribution originating from the term (t/texc)^ in (5.19) which 
asymptotically becomes tiny. As stated in the previous subsection, the actual Re(s^^'') spectrum is 
obtained by adding the above expressions over all excited KK-modes n. 
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The resonant peak rises at around k = kn, dividing wave modes by whether or not they went 
through the resonance band. Here, since we have adopted the simpUfying assumption of instan- 
taneous resonant amphfication, the above approximate expressions give a peak that rises suddenly 
at k = kn- Actually, the resonant oscillations arise sharply but within a finite wave number range, 
reflecting the finite time interval that a wave mode stays in the resonance band. As is shown in the 
first term in (C.6), when in the resonance band, s^^^ is parametrically amplified roughly proportional 
to (r — Tres). The rising region Ak around kn consists of wave modes for which the KK excitation at 

1/2 

Texc happens to be within the range Ar ~ — exc Tres (^) (5.16) around their resonance times Ties{k)- 
One can estimate that the resonant peak in the fc-space rises sharply within a wave number interval 
of 

AA: ~ eJLkn, (5.23) 

around k = kn- Thus actually the resonant oscillations peak at wave numbers slightly larger than kn- 
One can also see that in the rising region (5.23), effects from parametric resonance (i.e. the first line 
of (5.22)) is dominant over other effects (i.e. second line of (5.22)). 

In Figure 2 we illustrate a typical oscillation form induced on the curvature perturbation spectrum 
from an excited KK-mode with 9n ~ 0, i.e. Vnexc ~ 0. Wave modes that are sub- horizon at the 
KK excitation obtain oscillations about the offset — 2 fjne^dn exc cos 20 n- For modes in the range 
— T"exJ < k <kn which do not cross the resonance band, the oscillation period is a constant — vr/rexo 
and the oscillation amplitude is comparable to the offset. 

1/2 

The oscillation amplitude obtains a peak of ~ ??raexcCnexc at the wave number k ^ kn, be- 
yond which the amplitude decays proportionally to oc k~^, with A;-dependent oscillation period 
—ir/T^esik) ~ inexck- Here one clearly sees that the induced oscillations are determined basically by 
the two parameters f/^exc and ^nexc- 

Let us further remark that for nonzero On, oscillations with period — '^/'^ey.c is dominant also 
at wave numbers sufficiently larger than kn, as can be seen from the second line of (5.22). (Such 
oscillations at large k do exist even for On = 0, though they are suppressed. See discussions below 
(5.19).) On the other hand, oscillations in the range —t~J^^ < k <kn are suppressed for cos 20n ~ 1/4. 

For smaller ^nexo i-e. heavier KK modes, the resonant peak rises more sharply, and thus the 
envelope of the oscillations takes a more asymmetric form in fc-space. 

We have also numerically computed the Re(s^^'') spectrum as a function of the wave number, for 
the parameter set chosen in the previous subsection. The results are shown in Figure 3, where the 
blue solid lines show numerically calculated values of Re(s[.^^) at a time when all the displayed wave 
modes are well outside the horizon (and also when the KK mode has sufficiently damped away), 
and the red dashed lines denote the analytic estimations (5.20) - (5.22). One sees that the analytic 
result captures the overall behavior of the spectrum, though overestimates the resonant peak due 
to the approximation of the instantaneous parametric amplification. One can also see that as the 
dominant oscillation with period —Ti/T^cs{k) decays away for large k, tiny oscillations with period 
— tt/tcxc show up. Such oscillations are suppressed from Oi = 0, but are present at sub-leading orders 
beyond the approximate estimation of (5.22). 
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Figure 2: Schematic of the oscillations on the curvature perturbation spectrum induced by an excited 
KK-mode with On ~ 0, i.e. V'nexc ~ 0. The envelope of the resonant oscillations is sharper and more 
asymmetric for smaller ^nexo i-e. for heavier KK modes. 



Before ending this subsection, let us repeat that the particular forms of the oscillations generated 
by the KK modes can be attributed to a few parameters, which are r/„exc representing the excited 
KK amplitudes, and ^nexc setting the effective masses of the KK modes. In addition to these two 
parameters, the time of KK excitation Toxc determines where in A;-space the oscillations rise. 

5.4 Effects from Multiple KK Modes 

As we have stated above, for small KK excitations, effects from multiple KK modes simply add up. 
As an example, in Figure 4 we show the numerically computed Re(s[.^^) spectrum for the case where 
two KK modes n = 1, 2 are excited. The parameters are the same as in the previous subsections, 
except for that the n = 2 mode is further excited with a2 = 200, ■02 exc = 0.0004Mp, 62 = 0. (Note 
that an oc n, ^„exc oc n~^, f/„exc oc n^Vnexc) ^'^d kn oc n.) The plot shows superimposed oscillations 
from the two KK modes, which are peaked at around ki ~ A2kp and ^2 ^ 85A;p. 

Let us also note that effects from multiple KK modes can cancel each other, given that the 
excited KK amplitudes obey specific hierarchies such that different n modes produce Re(,s^^^) with 
similar oscillation amplitudes/offsets, and further if the phases On take appropriate values. When 
focusing on a certain wave number k and considering the time evolution of Re(sj3''), its oscillation 
amplitude can be suppressed (instead of amplified) when going through multiple resonance bands, 
if the phases On of the KK modes are suitably distributed, cf. first line of (5.19). However, the 
resonant oscillations of the Re(s^^^) spectrum in /c-space (i.e. first line of (5.22)) cannot be cancelled 
entirely, since ^nexc oc and thus different n modes source different oscillation periods. 

As we have seen explicitly, the KK tower generates oscillations on the curvature perturbation 
spectrum with specific patterns characterized by the integers n. It would be very interesting to inves- 
tigate such specific patterns in cosmological observables, which can allow us to extract information 
about the extra dimensional space. 
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Figure 3: Asymptotic values of Re(s^ ), denoting the KK induced oscillations on the curvature 
perturbation spectrum \{Vc_{k) /V(^Q{k) — 1), as a function of wave number in linear (left figure) and 
log (right) scales. Blue solid lines show numerically computed results, while the red dashed lines are 
the analytic estimations (5.20) - (5.22). 
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6 Strong Resonance from Large KK Excitations 



In this section we study cases with large KK excitations beyond the condition (3.20), i.e., when (3.7) - 
(3.10) are no longer good approximations of the inflationary dynamics. We will see that when the KK 
mode amplitudes are as large as ?/nexcCraexc ^0.1, the inflaton field fluctuations are strongly amplifled 
through parametric resonance, generating sharp spikes in the curvature perturbation spectrum. 



6.1 Curvature Perturbation Spectrum 

Let us first show how the resulting curvature perturbation spectrum changes when the KK mode is 
excited with larger amplitude. In Figures 5, 6, and 7, we show the results for the case where the 
n = 1 KK mode is excited at about 50 e-foldings before the end of infiation. The parameters are 
chosen similarly as in the previous section, i.e. (/)exc ~ 8.2Mp, fj, ~ 0.0016Mp, ai = 100, but with 
various KK mode amplitudes V'lexc- We again use kp to denote the wave number that exits the 
horizon at around the time of KK excitation. (The following parameters are also the same as in the 
previous section: eexc ~ 0.015, ^lexc ~ 0.082, and ki « 42A;p, though ki defined in (5.14) no longer 
denotes the position of the resonant peak for strong resonance.) The excited KK mode amplitudes 
are taken as ■0iexc = O.OOlMp (771 exc ~ 0.0018), 'tpiexc = O.OlMp (r?iexc ~ 0.18), and ipiexc ~ 0.033Mp 
(f?iexc ~ 2.0), respectively, for Figures 5, 6, and 7, with vanishing initial velocity tpiexc = 0. Upon 
obtaining the power spectrum, we have numerically solved the full set of equations of motion for the 
homogeneous background (3.2), (3.4), and (3.6), as well as for the infiaton field fiuctuation qk (4.5).^^ 
Note that even when the KK modes are largely excited, the power spectrum can be obtained by 
evaluating (4.43) when the wave modes have exited the horizon (4.41) and the KK mode oscillations 
are damped away such that r]n is sufficiently smaller than unity, satisfying for e.g. (4.42). This 
gives the power spectrum within errors of 0{e, r]n) upon evaluation. We have assumed the KK mode 
to be excited suddenly to tpnexc, which sources unphysical jumps at texc in quantities such as the 
Hubble parameter. However, we expect this treatment to be valid for discussing the KK signals 
since it is the parametric resonance following the KK excitation that sources sharp features in the 
perturbation spectrum. 

In the figures on the left side, the blue solid lines denote the power spectrum "P^ of the curvature 
perturbations, shown with log scale for the wave number. The black dot-dashed lines represent the 



^''If the initial conditions of qk were set by the Hankel type solution (4.24), it would contain errors of 0(e) (in the 
absence of KK-modes). Since we are dealing with a non-canonical inflaton, such errors modify the sound speed and 
induce artificial oscillations of 0{e) on the resulting perturbation spectrum, which can be confused with the KK effects. 
Therefore in our numerical computations we have set the initial conditions using the WKB-type solution: 

^'= = -(2^)3 (2^.(r))V2"H"'/ ^'^^'^'^^'} 

where 

™k{rj - 

This is a solution to the equation of motion (4.5) while the wave modes are well inside the horizon, within errors (in 
the absence of KK-modes) of 0{e^^, (^)'*), providing more accurate initial conditions than (4.24). One can also 
check that (6.1) satisfies (4.22) within errors of 0{e), and approaches qk oc e~'*'^ in the asymptotic past. (The 0{e) 
error for satisfying (4.22) can affect the overall amplitude of the perturbations, but does not directly lead to sourcing 
artificial oscillations in the spectrum.) 



(6.2) 
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spectrum "P^q in the absence of the KK mode. On the right side, we show in blue sohd lines the ratio 
("P^ — 'P^o)/^CO denoting the relative amplitude of the KK mode effects, with linear scale for the wave 
number. For comparison, we have also plotted the red dashed lines which are the asymptotic values 
of 2Re(s^^^), obtained by solving (4.33), with (4.25), (4.30), and (4.31). As we have discussed in the 
previous section, 2Re(s^^'') represents {V( — 'P(o)/'P(^o in the small KK regime fjnexc ^ 1 (cf. (4.46)), 
as is seen in Figure 5.^^ However, 2Re(s^^^) fails in describing the KK effects as fjnexc increases and 
the small KK approximations break down. As the KK amplitude increases, the resonant oscillations 
around k = ki grow asymmetrically towards the positive side of "P^Oj and eventually form sharp 
spikes on the power spectrum. 

The strong parametric resonance significantly enhances the inflaton field fluctuations, and the 
resulting spikes in the perturbation spectrum grows non- linearly with respect to fjnexc- A slight 
increase in the KK amplitude enhances the spike by orders of magnitude, for e.g., ij^nexc ~ 

0.05Mp 

sources (P^ — P(^o)/^fo as large as ~ 3000 at its peak. It should also be noted that the peak of the 
oscillations/spikes is shifted slightly towards larger k for strong parametric resonance, as can be seen 
in Figure 7. While the amplitude of the oscillations on P^ are significantly enhanced by the strong 
resonance, the oscillation period follows similar behavior as for the weak resonance. Moreover, the 
oscillations at wave numbers away from the peak (e.g. k < 40fcp and k > 70/cp in Figure 7) are not 
so different from the weak resonant calculations represented by 2Re(s^^''). Here we also remark that 
when the resonance is very strong, the spectrum can obtain multiple peaks, although the one at the 
wave number slightly larger than ki is the largest. 

The resonant spikes grow significantly as the KK amplitude is increased, until ipnexc becomes so 
large that the field dynamics are completely modified and/or the higher order KK interactions are 
important. For the Lagrangian (2.13) with the set of parameters chosen here, the KK mode no longer 
oscillates when tpnexc 

> 0.06 (or 

fjnexc ^6). In addition, i?^^ becomes negative and thus the inflaton 
field fluctuation q is tachyonic, cf. (4.3). However we should remark that in such case the product of 
fjnexc and (,nexc IS as large as fjnexdnexc ^ 0.5 and the condition (3.24) is almost broken. Hence the 
cubic or higher order KK terms in the Lagrangian can become important, and may prevent q from 
becoming tachyonic. It would be very interesting to explore this regime, where the higher order KK 
terms may contribute to source even stronger resonant features in the primordial spectrum. 

6.2 Effective Frequency 

The behavior of the strong resonance can be better understood by studying how the oscillation 
frequency of q^ resonates with that of the KK modes. For this purpose, let us redefine the field 
fluctuation as 

Uk = {aB^^Y^'^qk, (6.3) 
so that its equation of motion takes the form 

= Ufc + /|cff^fe> (6-4) 

^^The slight difference in the center of the oscillations between 2Re(s^^-') and {Vq — 'Pqo)/'Pqo arise due to the sub- 
leading KK contributions that were dropped in the s^^' computations. For example, the 0(7)„^„) corrections to the 
conformal time (4.11) slightly affects the expansion of the universe and thus modifies the offset of the perturbation 
spectrum. Such KK mode effects on the expansion history also source the deviation of Vc^ from Vi^o for wave modes 
that are super-horizon at the KK excitation, as can be seen in the figures on the left side. 
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where f^^g denotes the effective frequency squared defined as 



+ 




(6.5) 



The oscihating KK modes force the frequency /l^g itself to oscillate in time, and when becomes 
similar to the KK mode frequency ~ ttikk, the field fluctuation Uk experiences parametric resonance. 
The weak resonance discussed in Section 5 corresponds to tiny oscillations of f^^g, and the strong 
resonance happens when f^^g oscillates considerably. 

In Figure 8 we plot the time evolution of the field fluctuation and its effective frequency 
/fceff case of ipiexc ^ 0.033Mp and k ~ 57kp. This wave number is where the resonant 

spike is peaked at in Figure 7. The time is shown in units of the KK mode oscillation period 
27r/mKK — '^4'exc/ctiVcJc therefore is equivalent to the number of oscillations of V'l- Here, recall 
that the KK modes affect f^^g through quadratic terms (see e.g. (4.6)), hence f^^g oscillates with 
half the period of the KK mode. The KK mode is excited at t = in the figures. In the left figure 
we show the growth of the fluctuation \uk\ relative to its value at the KK excitation |nfcexc|- In the 
right figure, the blue solid line denotes the ratio /|eff/"^KK between the effective frequency of Uk 
and the KK mass. Moreover, the red dashed line shows k"^ /a^m^^, while the black dot-dashed line 
at unity represents the KK mass squared. We remark that the displayed wave mode k w 57kp is 
inside the horizon in the plotted time range. 

In the absence of the KK modes, the effective frequency /l^g traces k"^ /a? when the wave mode is 
inside the horizon. This is still the case for small KK excitations triggering weak resonance discussed 
in Section 5. However when the KK modes are largely excited, /Igg considerably oscillates and can 
even become negative. From the figures, one sees that Uk is enhanced by an order of magnitude 
soon after the KK mode excitation, while /|gff wildly oscillates around rn^-^. The fluctuation uj. 
initially oscillates with twice the frequency of the KK mode (i.e. the period matches with the tics on 
the time axis), and is enhanced at each oscillation. As /|gg becomes smaller than m^-^, parametric 
resonance ceases and \uk\ starts to increase slowly in time as oc~ a^/^. (Note that, in the absence 
of the KK modes, \uk\ scales proportional to a^/^ when inside the horizon, and as a?/"^ after horizon 
exit, as can be seen from its definition (6.3) and (4.24).) Let us also remark that the parametric 
resonance of the inflaton field fluctuation shares common features with preheating [38, 39]. We have 
studied that parametric resonance happens for Uk when its effective frequency /|gg oscillates with 
frequency similar to itself, hence in this sense, both the weak and strong resonance is analogous to 
the narrow resonance in preheating. 

In Figures 5 - 7 we have seen that the resonant peak in the perturbation spectrum slightly 
shift towards larger k for stronger resonance. This can be understood from the oscillating /|gff 
mostly taking values that are smaller than fc^/a^. In other words, the characteristic frequency of 
Uk effectively becomes smaller than k/a, and therefore the resonant features shift towards (slightly) 
larger wave numbers. 

The oscillations of the effective frequency /|gg are quickly suppressed as the KK mode oscillations 
damp away. The wave modes that enter the resonance band while the oscillations are large experience 

^^Now we are beyond the small KK approximation, but the KK mode ip\ still oscillates with period ~ 27r/mKK — 
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strong parametric resonance, and moreover, such resonant effects depend non-linearly on the KK 
amplitude (squared). This is why the strong enhancement of the curvature perturbations can be 
locahzed to a narrow wave number range, forming spikes in the perturbation spectrum. 




2 4 6 8 27r 



Figure 8: Time evolution of the inflaton field fluctuation Uf. (left figure) and its effective frequency 
squared (right figure), for ^lexc ~ 0.033Mp and k ~ 57kp. The time is shown in units of the 
oscillation period of the KK mode 27r/mKK5 ^-nd the KK mode is excited at t = 0. 



6.3 Condition for Strong Resonance 

For the set of parameters chosen in Section 6.1, the oscillation of the effective frequency /l^g (6.5) 
is dominated by its components 

and furthermore, by the same terms that were important for weak resonance, i.e., the terms ex- 
plicitly written in (4.6), (4.7), and (4.8). Recall that these terms originate from the three kinetic 
couplings (4.14) in the inflaton Lagrangian. Here we focus on those terms and quantitatively discuss 
the strong parametric resonance. We note that our aim here is to give a rough estimate on when 
strong resonance happens, and in this subsection we use some identities without proof. 

As a rule of thumb, we consider strong resonance to happen for wave modes in the resonance 
band if /|gg oscillates with amplitude as large as k'^/a^. This happens if, for e.g., B^p^ in the 
denominator of the first term in (6.6) oscillates substantially. Using ipn — 4>n cos{2anV^/'^t / (p + • • • ) 
and ipn — —2anV^^'^^pnSin{2anV^/'^t/(j) + •••), then one finds from (4.6) that i?^^ oscillates with 
amplitude of order unity or larger for (for simplicity we consider a single KK mode excitation) , 

0.l<al^=fine.c^ne.c. (6.7) 

Vexc 

One can further check that the other terms in (6.6) also oscillate /|gg with amplitude ~ fc^/a^ for 
similar KK amplitudes. Therefore, we can consider (6.7) as the rough indicator for strong resonance 
to happen. 

A wave mode undergoes strong resonance if it approaches the resonance band (i.e. k? ja^ ~ "'-kk) 
while (6.7) holds. Denoting the largest (smallest) wave number that experiences strong resonance 
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by k 

max (^min); then from. fin(,n oc~ 0. One can estimate the ratio 

~ (lOf/nexC^nexc)^^^ • (6.8) 

"•min 

Therefore, in the regime (3.24) where one can neglect the cubic or higher order KK terms in the 
Lagrangian, the wave interval over which strongly resonant features arise is quite narrow. Moreover, 
since the strong resonant amplification non-linearly depends on the KK amplitudes, the resonant 
features can further have hierarchical structures. 



7 Conclusions 

Wrapped brane infiation necessarily possesses infinite KK degrees of freedom in addition to the zero- 
mode inflaton, and the brane's Nambu-Goto action naturally gives kinetic couplings among them. 
In this paper, we have investigated resonant amplification of the curvature perturbations triggered 
by the KK oscillations during inflation. We found that the resonant signals in the perturbation 
spectrum can be localized to narrow wave number ranges, in contrast to previous studies where 
signals were normally distributed over a rather wide A:-range. In the effective four dimensional 
theory, the zero-mode inflaton couples to the heavy KK modes which are excited and oscillates 
during inflation. Small KK excitations lead to weak resonance, sourcing oscillatory features in the 
perturbation spectrum (cf. Figure 2, 5) whose oscillation amplitude is proportional to the excited 
KK amplitude squared. For larger KK excitations, the resonant amplification becomes extremely 
efficient with non-linear dependence on the KK amplitude. We saw that such strong resonant effects 
damp more rapidly than the KK oscillations during inffation, therefore leave spiky features localized 
to a narrow wave number range in the perturbation spectrum (cf. Figure 7). Both weak and strong 
resonant effects were mainly sourced by the couplings between the inflaton kinetic terms and the 
oscillating KK modes. 

Motivated by the string theory construction of [6], we have supposed the brane tension to be 
the main driving force for inflation and considered only the Nambu-Goto action. However we note 
that, even if there are further contributions to the brane potential, the Nambu-Goto action gives the 
kinetic terms that contain the important interactions with the KK modes. Therefore we expect the 
resonant effects studied in this paper to be rather generic features for inflation driven by extended 
sources. From an effective theory point of view, it would also be interesting to systematically examine 
resonant effects arising from kinetic couplings in general. 

We have restricted ourselves to KK modes that are not excited too largely and thus studied up 
to quadratic KK interactions in the Lagrangian. Some region of the strong resonance was close to 
the regime where the higher order KK terms can become important, thus it is of interest to analyse 
effects from higher order interactions. Especially when going beyond the strong resonance regime 
we have studied, the tower of KK interactions may further enhance the resonant signals. 

The wrapped brane model provides an example where the kinetic couplings between the inflaton 
and the heavy KK modes produce sharp features in the curvature perturbation spectrum, and to 
that end we have focused on the evolution of the perturbations under the KK oscillations. The 
KK modes can be excited as the wrapped inflaton brane passes by structures or localized sources 
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along the inflationary trajectory,^ however we have not yet specified the excitation mechanism in 
detail. It would be very interesting to study the whole picture including dynamical excitation of the 
KK modes. One can imagine that the KK excitation is triggered as the inflaton comes across some 
point(s) in field space. Then the KK excitation itself may become inhomogeneous and give further 
contributions to the curvature perturbations. It could also be that the KK modes are repeatedly 
excited during inflation, if, for e.g., the trajectory of the wrapped brane is a repeated circuit as in 
monodromy models [6, 7]. In such case, the spikes/oscillations would repeatedly show up in the 
perturbation spectrum with the interval corresponding to the size of the inflaton circuit. We also 
note that, when constructing explicit models for wrapped brane inflation, the fate of the excited KK 
modes need to be taken into account, especially in relation to the moduli problem [41, 42, 43]. (See 
also [44] for discussions in this direction.) 

Interestingly, the resonant signals in the curvature perturbations are characterized by the mi- 
crophysics of the KK modes, which is tied to the properties of the internal manifold. The height of 
the spikes/oscillations are basically determined by the excited amplitude of the KK modes, while 
the width and oscillation period of the features in k-space are set by the KK mass corresponding to 
the size of the wrapped cycle. The positions of the resonant signals in /j-space tell us when during 
inflation the KK modes were excited (possibly denoting the positions of localized sources in the 
internal manifold), and furthermore, the KK tower generates signals periodic in wave number k, cf. 
Figure 4. All these features open up possibilities for probing the extra dimensional space through 
cosmological experiments. It will also be interesting to study how well we may detect and measure 
these resonant signals through observations of the cosmic microwave background (CMB) and large 
scale structure. The oscillatory features from weak resonance may be somewhat smoothed out by 
the CMB transfer function (as is discussed in, e.g. [25]), but the spikes from strong resonance are 
expected to leave distinct imprints in the CMB temperature anisotropy. We leave to future work a 
detailed analysis of the observational consequences of weak/strong resonance during inflation. We 
should also mention that the resonant signals can show up in the non-Gaussian signals as well, as 
discussed in [16, 17, 23]. It is quite possible that the strong resonance which sources spikes in the 
power spectrum also leave violent marks on the higher order correlation functions. 
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A Primordial Density Perturbations from Multi-Field Inflation with 
Various Kinetic Terms 

In this appendix we consider density perturbations from multi-field inflation with an action contain- 
ing various forms of kinetic terms. We make use of the (JA/'-formalism [33, 34, 35, 36] for obtaining the 

^^See also [40] which discusses KK excitations due to oscillations in the compactification scale of the extra dimension. 
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density perturbations, so here we simply derive the action of the field fluctuations up to quadratic 
order. The calculations are an extension of the previous works discussing inflationary models with 
multi-fields and/or non-canonical kinetic terms, see for e.g. [45, 46, 47, 48, 49, 50, 51]. 
The action we consider is of the following form: 

S = I d^x^g I + L (/, X^) I , (A.l) 

where ip^ denotes the fields (labeled by /), and the kinetic terms, 

X^ = -Gfjg^^'^d^p'd,^'. (A.2) 

Gfj is the field space metric, where A labels different metrics, i.e. different forms of kinetic terms. 
We impose Gfj = Gjj, and consider Gfj to be constants (in other words, independent of (f^). Note 
that in this appendix we label fields by /, J, • • • , and field metrics hy A, B , ■ ■ ■ . Moreover, we express 
partial derivatives in terms of ip^ and X"^ as, respectively, 

dL ^ dL , 

Then one can derive the energy-momentum tensor, 

dL 

T^u = 9f.uL - 2— = g^,L + 2LAGfjd^ip^ d^if^ (A.4) 
as well as the equation of motion of ip^ , 

Li + -^d^ [V^LAGfjg^'^d^ip') = 0. (A.5) 

A.l Homogeneous Background 

First let us derive the equations of motion of the homogeneous background. We fix the background 
metric to a flat FRW: 

ds^ = -dt^ + a^{t)5ijdx'dx^ , (A.6) 

where i, j = 1, 2, 3 run over the spatial directions. 
Then the Einstein equations are 

3M2//2 = -L + 2LAGfj^^ip\ (A.7) 

- M^{2H + 3H'^) = L, (A.8) 

where an overdot denotes a time derivative, and H = d/a. 
The equation of motion of 93^ (A.5) now takes the form 



Lj - 2 [LAGfjif') ' - GHLAGfjif' = 0. {AS 
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A. 2 Field Fluctuations 

We then study fluctuations around the homogeneous background. We adopt the ADM formahsm, 

= -N'^dt^ + hij{dx' + N'dt){dx^ + NHt), (A.IO) 
under which the action (A.l) is rewritten as (we take > 0): 




(A.ll) 



Here h = det(/iij), and R^'^^ is the scalar curvature on the t = const, spatial hypersurface. The 
symmetric tensor Eij is defined as 

Eij = \ (hij - ViNj - VjNi^ , (A.12) 

where Vj is a derivative associated with hij, and E = h^^Eij. We also note that the indices i,j are 
raised and lowered by /i*-' and hij, respectively. is now expressed as 

r 1 / /V* JV^ \ /V* 1 



Hereafter we take uniform curvature slicing, such that t = const, slices have vanishing Ricci 
curvature, thus 

hij = a{tf5ij, =0, Vi ^ di, (A.14) 



and the action becomes 



S = / dtd^xVh 



(E,iE'^ - E^) +NL 
2N ^ ^ ' 



We will focus on the field fluctuations defined as 



(A.15) 



(A.16) 



where ip^ is the homogeneous classical background, and later on convert them into the curvature 
perturbations using the (^A/"- formalism (A. 36). 

The lapse N and shift A* are Lagrange multipliers in the action (A.15), hence their equations of 
motion can be used as constraints, which are, respectively, 



N^L = {E.jE'i - E^) + 2LAGf,yv\ 



LAh'^GfjV^djif^ 



N 



Here we have defined 



= - A*5,/. 



(A.17) 
(A.18) 

(A.19) 
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Let us rewrite the lapse and shift as 

N = l + a, 

(A.20) 

Ni = dii^ + Ni, 

where Ni is the incompressible part, i.e., 

h'^djNi = 0. (A.21) 

We further expand the functions in terms of the field fluctuations Q, 

a = ai + a2 + ■ ■ ■ , 

= V'l + ^^2 + • • • , (A.22) 
N, = N^^^ + N^^^ + -- - , 

where the numbers represent orders of Q, e.g. On = 0{Q'^). We require the constraint equations to 
hold order by order. 

In order to expand the action up to second order in Q, one needs to consider the lapse and shift 
only up to the first order, i.e. qi, ipi, and N^^^ (cf. [45, 48]). The constraint equation (A. 17) at its 
zeroth order reproduce the Friedmann equation (A. 7), while its first order is 

H 



- AL^X'GfjQ'^i - 2 [L^XiQ' + LZ^r) Gfj^i^i, (A.23) 

where we have used (A. 7) upon obtaining this form. The zeroth order of (A. 18) is trivial, and the 
first order is ^ 

2d, {ulHai - LfGfj^iQ') = ^h^>^djdkN^'\ (A.24) 
bmce 0di{h^^dkdiNf^) = 0, after choosing proper boundary conditions, one arrives at 

h'^didjN]^'^ = 0, (A.25) 

MlHai = LfGfj^iQJ. (A.26) 

We also note that the evolution equation (A. 8) is obtained by combining the zeroth order of (A. 17) 
(i.e. the Friedmann equation) and the zeroth order equation of motion of (f^. 

Using the zeroth order of (A. 17), the zeroth order equation of motion of (p^ , and the first order 
of (A. 18) (i.e. (A.25) and (A.26)), one can expand the action up to second order in Q, 

S = j d^xdta^ -3M^H^ + L(°) - if Cfj^d'^ diQ^ djQ-^ 

+ {^^5°^' + il^^ {$LfLfGf,Gj,^^^^^]Q^Q^ (A.27) 
+ [if Of J + 2L%Gfj,G^j^^^^^,] Q'Q' + 24°]g^^^^Q^Q^ + 0{Qf 
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where "=" is used to denote that we have dropped total derivatives in the integrand. 

Hereafter we omit the sub(super)script denoting the homogeneous background. Further intro- 
ducing 

p 

Bij = 2LAGfj + ALABGfKG^L^'^ip 
Cij = 2LaiG^j^^ , 

(note that Mjj = Mjj and Bjj = Bjj , but C/ j is not necessarily symmetric) , then the second order 
action for Q can be rewritten in the following form: 



Mi J = Li J + ( -LALnGf^G^Lip^ ^p' 



A r^B ,.K,.L (A.28) 



Si = / dtd^xa^ 



^LAGfjdQ'dQ^ + Umij-^ (a^C/j) ) Q^Q^ 



+ IbijQ'Q-' + ^{Gij - Gji)Q'Q' 



(A.29) 



where we use the abbreviated expressions 



dXdY = 6'^diXdjY, d'^X = 5'^d^djX. (A.30) 



The equation of motion for can be obtained from (A.29) as 



= ^LAGjjd^Q'^ + MjjQJ + 2LAiGjKip^Q'' - ^ [a^ [BjiQ-^ + CjiQ^) }' . (A.31) 

IS of the conformal time 

dt = adT. (A.32) 



The results can also be expressed in terms of the conformal time 



Redefining the field as 

q^ = aQ^, (A.33) 

and further Fourier expanding , 

q'{r,x) = d'fce-^-gi(T), (A.34) 

then the equation of motion (A.31) can be rewritten as 



= Bjiq'i^ +{Cji- Gij + Bji) aq'k^ 



+ ^2LAGjik^ - Bji^ + (-Mj, + Cji^j + (^-Bji + 2Cji + Cij^j a'j ql (A.35) 
Here k = \k\, and a prime denotes a derivative with respect to the conformal time r. 

These are the main results of this appendix, which can be used to obtain the field fluctuations 
after choosing appropriate initial conditions (which corresponds to choosing the vacuum state) . Then 
one can calculate the resulting curvature perturbations using the (5AA- formalism [33, 34, 35, 36], 

C = AfjQ' + ^AfijQ'Q^ + ■■■ , (A.36) 
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where J\f is the number of e-folds between an initial flat hypersurface and a final uniform density 
hypersurface. The right hand side of (A. 36) can be computed at any time after the separate uni- 
verse picture becomes a good approximation, provided that there are no isocurvature perturbations 
sourcing further SAf. The computations of the curvature perturbations using (A. 35) are carried out 
in detail in Section 4, focusing on the specific action of (2.13). 



B Homogeneous Functions in the Field Fluctuation Action 

Here we write down the forms of the homogeneous functions defined in Appendix A (at e.g. (A.28)), 
when applied to the action (2.13). They show up in the second order action of the field fluctua- 
tion (4.3). 

Expressing the kinetic terms as 

X = -{dcl)f, X = -Y^{d^Pnf, Xn = -{d<t)-d^n) (forn ^ 0), (B.l) 
then the Lagrangian (2.13) can be written as £ = ^—gL where 



L 



+ ^yfx^ + i4^2-^^XX„V (B.2) 



with 



-1/2 

7=(1-^ ■ (B.3) 



With the above L, the functions in the fluctuation equations are expressed as follows: 

LAGi^ = ^. (B.4) 

= 2^ + A—^4>' + 4 X: Q^Hn + E QjlJ^Jm^n, (B.5) 

''^^ = 'W^^^T.^Jn. (B.7) 
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C Computation of Weak Resonance from Oscillatory Sources 



In this appendix we give detailed computations of weak resonance from small KK excitations. We 
focus on the equation (5.6) and study the inflaton field fluctuations induced by oscillatory source 
terms. Solutions of the evolution equation of the form (5.10): 

- + '^ikspir) = i^sin(a;r + B) (C.l) 

are analyzed, where k, oj, and are real (here k and uj are positive), while A is a complex number. 
We start our discussions by supposing that these parameters are constants. We arbitrarily set the 
initial conditions at a certain time 



The homogeneous solution of (C.l) is (5.2), and the particular solutions are given in (5.11) and 
(5.12). 



C.l Solution for 2k ^ to 

When 2k ^ uj, the solution of (C.l) with initial conditions (C.2) is 



4'V) 



A 



2k 



CO 



cos{uJT + 0) + 2 sin(a;r + 0) 



2k 

OJ 



UJ 

2k 



cos(a;r* + G) 



\^ism{ojn + 0) + ^ cos(a;r, + 9)} e^^'^^— *) 



+ ^* + 2k' 



/ _g2ifc(r-r.)| (C_3) 



Especially when s* 



and further 2k S> uj, then the solution is 



20 



A 

'2kuj 



UJ 



cos{ujT + 0) + i— sin(ci;T + 0) — cos(a;r* + 0) 
2k 



i|^sin(a;r*+0)e2*'=("-"*)+O 
Ik 



UJ 



while s* 



and 2k <^ UJ give 



4"M 



A 
2kuj 



cos wr. 



.^ + e){i-e2^^(---)} + o(^) 



(C.4) 



(C.5) 



One sees from (C.4) and (C.5) that the system basically chooses the smaller of the two frequencies 
k/iT and uj/2tt. This can also be understood as follows: In the case 2k ^ uj, the particular solution 
(5.11) mainly determines the oscillatory behavior, while for 2k <^ uj the homogeneous solution (5.2) 
dominates. Therefore, when further taking into account the time variation of the parameters such 
as A, then (C.4) shows oscillations with a time-varying amplitude, while (C.5) gives a constant 
oscillation amplitude. 

However, it should also be noted that the story is more complicated when focusing only on the 
real or imaginary part of s^^K Especially for the 2k ^ uj case (C.4), the amplitude of oscillations 
with frequency a;/27r for Re(s[.^^) changes drastically depending on whether A is real or imaginary. 



''Here we use x — 0{y) for a complex variable x to denote < \y\. 
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C.2 Solution for 2k = uj 



For the resonant case 2k = lo, the solution for arbitrary and s'^ is 

2ik{T - T,)e^(2^"+®) - cos(2A:r + 6) + {2 - e^^^^"""*)} cos(2A;r, + 9) 

+ .. + ^.:{l-e2^'=(--)}. (C.6) 

The term 2ik{T — r*)e*^^'^'^'''®'' in the first hne represents parametric resonance, and since we have in 
mind wave modes weU inside the horizon (i.e. —kr ^> 1), one naively expects the term to dominate 
over the others in the [ ] parentheses. This can be verified as follows: Going back to the original 
equation shown in the first and second lines of (5.6), one should recall that parametric resonance 
happens while the characteristic oscillation with phase 2kT is synchronized with that of the KK- 
induced source term with a;*r* ln(T/T*) (where we have explicitly written c<j* to indicate uj at r*). 
Fixing the time t^, to be when 2k = uj^,, then one can estimate the time scale for the resonant 
solution (C.6) to be valid by calculating how long 2k{T — t^) « 2fcn ln(T/r^,) lasts. Considering the 
wave mode k to exit the resonance band when the two phases are misaligned by half-period vr, one 
can compute this time r/ as 

Tf-n=nln^-^. (C.7) 

Since —kr^ » 1 one obtains 

k{Tf - ^ {-7rkn)^/\ (C.8) 

which shows that the resonant term 2ik{T — r*)e*^^'^'^+®) does become much larger than the others 
in the [ ] parentheses of (C.6) by the time the wave mode leaves the resonance band. 
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C.3 Time Dependent Parameters 

We now extend the above discussions and consider the parameters A, oj, and in (C.l) to evolve in 
time, as in the original equation (5.1) from oscillating KK modes. The parameters are considered to 
vary with time scales much larger than the oscillation periods ir/k and 2tx jio. Moreover, we suppose 
io to monotonically increase in time r, and A to decrease. Here we provide rough expressions that 
track the time evolution of ^ from the KK excitation when we set the initial conditions 

4'^(rexc)=0, /i'^(rexc) =0, (C.9) 

until some time r (< — 1/A;). In other words, we focus on wave modes that are well inside the horizon 
upon the KK excitations, and study their sub-horizon evolutions. For the wave mode that is in the 
resonance band at the KK excitations, we refer to its wave number as i.e., 2kn = w(tcxc)- 

Keeping in mind the discussions in Section C.l, we obtain rough approximations for s^^\t) by 
collecting the leading contributions from the above solutions, under the following procedure: For 
the oscillations with frequency k/ir, we treat them as the homogeneous solution (5.2) and consider 
their oscillation amplitudes as frozen at the initial value. The offset of the oscillations are treated 
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similarly. Thus for wave modes — I/texc ^ k < kn that do not cross the resonance band, 2k < uj 
holds throughout and one sees from (C.5) that the time evolution of s^^^ is approximated by 

(1), N ^exc ^ „ ^ \^ „ _ 1 



COs(a;excTexc + eexc){l-e2*^("-"-)} for T,^c<T<--. (C.IO) 



Here, note that no matter Aexc is real or imaginary, the oscillatory term e2fcj(T-rcxc) gives similar 
oscillation amplitudes to both the real and imaginary parts of s^^\ Therefore we have neglected the 
0{k/uj) terms in the parentheses in (C.5), which are sub-leading anyway. 

On the other hand, oscillations with frequency uj/2tt arise from the particular solutions (5.11) 
and (5.12), with time-dependent oscillation amplitude and frequency. Thus for wave modes that ex- 
perience parametric resonance, i.e. k > kn, the approximate expression before crossing the resonance 
band (i.e. 2k > uj) is obtained from (C.4) as 

A 



+ COs(WexcTcxc + ©exc) + sin(a;cxcTcxc + 0exc)e^*^^'' 

for Texc < r < Tres, (Cll) 

where we define Tres as when 2k = a;(rres)- (We also denote parameters at Tres with the sub- 
script "res".) Focusing on Re(s^^^), it should be noted that for real A the oscillation with fre- 
quency uj/2'ir dominates over that with k/ir, but for imaginary A the oscillations with different 
frequencies can have comparable amplitudes. 

Upon obtaining the expression for s^^'' after crossing the resonance band (i.e. 2k < co), we 
succeed the offset and the oscillatory term with frequency k/ir from the "homogeneous solution 
part" in (C.ll) (i.e. the second line), but replace the "particular solution part" (i.e. the first line) 
by the resonant term in (C.6) combined with (C.8). Here for simplicity we treat the parametric 
amplification to happen instantaneously. Since the amplified term oscillates with frequency k/ir 
upon leaving the resonance band, from then it can be considered as a homogeneous solution with 
constant amplitude. We ignore the terms in (C.5) which may arise after leaving the resonance band, 
since their amplitudes are smaller than the oscillation amplitude of the resonant term. Therefore we 
obtain 

4'^M--^(--^--)^/^e^(^'=-+®-=) 

+ ^^J^ COs(cJcxcTcxc + ©exc) + sin(a;cxcTexc + ©exc)e^*''^'^~'^'=""^ 

for rres<r<-i. (C.12) 
k 

Here we have chosen the phase in the first line of the right hand side such that it smoothly connects 
with that of the first line of (C.ll), when ignoring the time variation of the parameters such as cj. 
(However, we should also remark that this choice of phase is just a simplified procedure taken to 
smoothly connect the oscillation phase through the assumed instantaneous parametric amplification.) 

Thus we have obtained rough approximations for s^^'' (C.IO), (C.ll), and (C.12). It is shown in 
Section 5 that these expressions capture the overall behavior of s^^^ when inside the horizon. 
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